Search for two-scale localization in disordered wires in a magnetic field 
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The supersymmetry technique [A. V. Kolesnikov and K. B. Efetov, Phys. Rev. Lett. 83, 3689 (1999)] 
predicts a two-scale behavior of wavefunction decay in disordered wires in the crossover regime from 
preserved to broken time-reversal symmetry. We have tested this prediction by a transmission 
approach, relying on the Borland conjecture that relates the decay length of the transmittance to 
the decay length of the wavefunctions. Our numerical simulations show no indication of two-scale 
behavior. 

PACS numbers: 72.15.Rn, 05.60. Gg, 73.20.Fz 
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In a remarkable paper fil , Kolesnikov and Efetov have 
predicted that the decay of wavefunctions in disordered 
wires is characterized by two localization lengths, if time- 
reversal symmetry is partially broken by a weak mag- 
netic field. Using the supersymmetry technique M it 
was demonstrated that the far tail of the wavefunctions 
decays with the length £2 characteristic for completely 
broken time-reversal symmetry — even if the flux through 
a localized area is much smaller than a flux quantum. At 
shorter distances the decay length is £1 = -kfa. It was sus- 
pected that previous studies by Pichard et al. || found 
single-scale decay because of the misguiding theoretical 
expectation of such behavior. This expectation was also 
the basis for the interpretation of the experiments by 
Khavin, Gershenson, and Bogdanov g on submicron- 
wide wires. 

The prediction of Kolesnikov and Efetov calls for a 
test by means of a dedicated experiment or computer 
simulation. It is the purpose of this work to provide 
the latter. We target the key feature of the two-scale 
localization phenomenon, which is the doubling of the 
asymptotic decay length at infinitesimally weak magnetic 
fields. 

Our numerical simulations are based on a trans- 
mission approach. We rely on the Borland conjec- 
ture H (believed to be true generally ||), that relates 
the asymptotic decay of the transmittance T with in- 
creasing wire length L to the asymptotic decay of the 
wavefunction ip(L). According to the Borland conjec- 
ture, the Lyapunov exponent a — — limL^oo \L~ 1 lnT 
is identical to the inverse localization length £ _1 = 
— lim£^oo L^ 1 In |?/>(L)|. Moreover, £ and a are self- 
averaging, meaning that the statistical fluctuations be- 
come smaller and smaller as L — ► 00. Our numeri- 
cal simulations show that the crossover from £ = £1 
to £ = £ 2 does not occur until the flux $^ through a 
wire segment of length £1 is of the order of a flux quan- 
tum $0 = h/e. For our longest wires (L > 150£i) the 
crossover according to Ref. (jj should have occurred at 
$^/$o — exp(— £/8£i) ~ 10 . We consider various pos- 
sible reasons for the disagreement with Ref. ["J (finite 
number of modes, anomalously localized states), but be- 



e 



-200 




-300 - 



-400 



100 

FIG. 1. Average logarithmic transmittance (InT) as a 
function of wire length L for the Anderson model with N = 5 
propagating modes. The two dashed lines have the slopes pre- 
dicted for preserved (/3 = 1) and broken (/3 = 2) time-reversal 
symmetry. From bottom to top the data corresponds to fluxes 
$ 5 /$ =0, 0.0005, 0.005, 0.05 (four indistinguishable solid 
curves), 0.5, 1, 2.5, 5, 10, 15, 20, 25, 40, 50, 75, 125 (two 
indistinguishable solid curves). The inset shows InT for an 
individual realization with $^ = |$o (solid curve) and the 
slope of the ensemble-averaged result (dashed line). 



lieve that none of these provides a satisfactory explana- 
tion. 

Our first set of results is obtained from the numerical 
calculation (by the technique of recursive Green functions 

f) of the transmission matrix t for a two-dimensional 
nderson Hamiltonian with on-site disorder. In units 
of the lattice constant a = 1, the width of the wire is 
W = 13 and the wavelength of the electrons is A = 5.1, 
resulting in N = 5 propagating modes through the wire. 
The localization lengths £1 = (N + 1)1 and £ 2 = 2NI 
are determined by the scaling parameter / of quasi-one 
dimensional localization theory, which differs from the 
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FIG. 2. Asymptotic decay length (solid circles) and 
weak-localization correction ST (open circles) as a function 
of flux for the N = 5 Anderson model. 




transport mean-free path by a coefficient of order unity 
H . The average of the transmittance T — tr t$ in the 
metallic regime, fitted to (T) = N(l + L/l)^ 1 , yields 
I = 65. This gives a localization length £i = 390 for 
preserved time-reversal symmetry (symmetry index (3 — 
1) and a localization length £2 = 650 for broken time- 
reversal symmetry {(3 — 2). 

Fig. [I] shows the ensemble-averaged logarithm of the 
transmittance (InT) as a function of wire length L for 
various values of the magnetic field B (or flux <J>£ = 
W£iB). We find a smooth transition between the the- 
oretical expectations for preserved and broken time- 
reversal symmetry. Most importantly, we find an asymp- 
totic slope s(B) = limL^oo L _1 {lnT) that interpolates 
smoothly between the values s — — 2/£i for B = and 
s = — 2/£ 2 for large B. There is no indication of a 
crossover to the slope s = —2/^2 for smaller values of 
B, even for very long wires (L > 150£i). According to 
the theory of Ref. Q], the crossover should occur at a 
length L cross given by 



(1) 



icross/6 = 81n(Vl2$ /47T$ ? ) + 0(1), 



which is well within the range of our simulations (£ C ross — 
14£i for $£ ~ 0.05<&o). The absence of two-scale behav- 
ior in the transmittance of an individual, arbitrarily cho- 
sen realization is demonstrated in the inset of Fig. fil, for 
$£ = ^$o- The self- averaging property of the Lyapunov 
exponent is evident. 

The asymptotic decay length £(B) — —2/s(B) is plot- 
ted versus magnetic field in Fig. @, together with the 
weak-localization correction 6T — T(B = 00) — T(B) 
at L = £1. For both quantities, breaking of time-reversal 
symmetry sets in when $^ is comparable to $0. The tran- 
sition from /3 = 1 to /3 = 2 is completed for $£ « lOO^o- 

Our second set of results is obtained from a computa- 
tionally more efficient model of a disordered wire, consist- 
ing of a chain of chaotic cavities (or quantum dots) with 
two leads attached on each side. This so-called 'domino' 
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FIG. 3. Same quantities as in Figs, hi and H, but now for 
the N = 50 domino model. In the upper panel, the magnetic 
flux parameter 5 = 0, 0.0001, 0.0002, 0.0005, 0.001, 0.002, 
0.005, 0.01, 0.02, 0.05, and 0.1. In the inset, 5 = 0, 0.00001, 
and 0.0001 (indistinguishable curves). 



model [9j is similar to Efetov's model of a granulated 
metal p ] and to the Iida-Weidenmiiller-Zuk model of con- 
nected slices |1CJ • The length L is now measured in units 
of cavities, and the mean free path 1 = 1. The scatter- 
ing matrices of each cavity are randomly drawn from an 
ensemble (proposed by Zyczkowski and Kus [[llj) that 
interpolates (by means of a parameter 5) between the 
circular orthogonal (j3 = 1, S = 0) and unitary ((3 = 2, 
(5 = 1) ensembles of random-matrix theory. The relation- 
ship between 5 and ${/$o is linear for i5<1. 

We increased the number of propagating modes to 
N = 50, because it is conceivable that the two-scale lo- 
calization becomes manifest only in the large N- limit, or 
that only in this limit the critical flux <!>£ for the transi- 
tion from £1 to £2 becomes <C $o- (I n the experiments of 
Ref. B N rs 10, so our simulations are in the experimen- 
tally relevant range of N.) Because of the much larger 
value of N, we restricted ourselves for larger values of the 
magnetic flux to L ~ 25£i, which should be sufficient to 
observe the localization length £2 for Q^/Qq > 10 -2 . For 
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FIG. 4. Logarithm of the average transmittance ln(T) as 
a function of wire length L for the N = 5 Anderson model 
at various values of the magnetic field (solid curves; from 
bottom to top, <E>£/"I>o = 0, 5, 25, 50, 125). The dashed curves 
are the theoretical prediction of Refs. [ |13| , fl4[ for zero and large 
magnetic field. 



smaller values of the flux, we increased the wire length 
to L ~ 100£i. The data is presented in Fig. 0. It is qual- 
itatively similar to the results for the N — 5 Anderson 
model. Instead of two-scale behavior, we only see a sin- 
gle decay length which crosses over smoothly from £i to 
£2 with increasing S. Again, the crossover of £ coincides 
with the crossover of the weak-localization correction, so 
there is no anomalously small crossover flux for the lo- 
calization length. 

The logarithmic average (lnT) is the experimentally 
relevant quantity since it is representative for a single 
realization (see Fig. [I], inset). The average transmit- 
tance (T) itself is not representative, because it is domi- 
nated by rare occurrences of anomalously localized states 
]12J . Since Kolesnikov and Efetov B studied the average 
of wavefunctions themselves, rather than the average of 
logarithms of wavefunctions, it is conceivable that their 
findings are the result of such rare occurrences. For com- 
pletely broken or fully preserved time-reversal symmetry 
the average transmittance is given by |13| 



ln(T) = -L/2& - § In Lfo + 0(1). 



(2) 



The order 1 terms are also known p3J14 and contribute 
significantly for L < 30£i. (This is the numerically acces- 
sible range, because anomalously localized states become 
exponentially rare with increasing wire length.) We have 
plotted the full expressions in Fig. H (dashed curves), 
together with the numerical data for the N — 5 Ander- 
son model. Again we find a smooth crossover between 
preserved and broken time-reversal symmetry. There is 
no transition with increasing wire length to a behavior 
indicative of completely broken time-reversal symmetry, 
even though the flux <&£ is much larger than required 
[according to Eq. (|l|)] to observe this crossover for the 
wavefunctions. 



In conclusion, we have presented a numerical search 
for the two-scale localization phenomenon predicted by 
Kolesnikov and Efetov B, with negative result: The 
asymptotic decay length of the transmittance is found 
to be given by £1 and not by £2, as long as the flux 
through a localization area is small compared to the flux 
quantum. How can one reconcile this numerical find- 
ing with the result of the supersymmetry theory? We 
give three possibilities. (1) One might abandon the Bor- 
land conjecture and permit the asymptotic decay length 
of the transmittance (Lyapunov exponent) to differ from 
the asymptotic decay length of the wavefunction (local- 
ization length). Since the Borland conjecture has been 
the cornerstone of localization theory for more than three 
decades, this seems a too drastic solution. (2) One could 
argue that the wires in the simulation are too narrow 
or too short — although they are in the experimentally 
relevant range of N and L. (3) One could attribute the 
two-scale localization phenomenon to anomalously local- 
ized states, that are irrelevant for a typical wire. We can 
not fully exclude these two remaining possibilities, but 
they would severely diminish the experimental relevance 
of the phenomenon. 

A discussion with P. G. Silvestrov motivated us to look 
into this problem. We acknowledge helpful correspon- 
dence with A. V. Kolesnikov and support by the Dutch 
Science Foundation NWO/FOM. 
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